\(\int \frac {1}{(d f+e f x)^2 (a+b (d+e x)^2+c (d+e x)^4)^2} \, dx\) [651]

Optimal result

Integrand size = 33, antiderivative size = 360 \[ \int \frac {1}{(d f+e f x)^2 \left (a+b (d+e x)^2+c (d+e x)^4\right )^2} \, dx=-\frac {3 b^2-10 a c}{2 a^2 \left (b^2-4 a c\right ) e f^2 (d+e x)}+\frac {b^2-2 a c+b c (d+e x)^2}{2 a \left (b^2-4 a c\right ) e f^2 (d+e x) \left (a+b (d+e x)^2+c (d+e x)^4\right )}-\frac {\sqrt {c} \left (3 b^3-16 a b c+\left (3 b^2-10 a c\right ) \sqrt {b^2-4 a c}\right ) \arctan \left (\frac {\sqrt {2} \sqrt {c} (d+e x)}{\sqrt {b-\sqrt {b^2-4 a c}}}\right )}{2 \sqrt {2} a^2 \left (b^2-4 a c\right )^{3/2} \sqrt {b-\sqrt {b^2-4 a c}} e f^2}+\frac {\sqrt {c} \left (3 b^3-16 a b c-\left (3 b^2-10 a c\right ) \sqrt {b^2-4 a c}\right ) \arctan \left (\frac {\sqrt {2} \sqrt {c} (d+e x)}{\sqrt {b+\sqrt {b^2-4 a c}}}\right )}{2 \sqrt {2} a^2 \left (b^2-4 a c\right )^{3/2} \sqrt {b+\sqrt {b^2-4 a c}} e f^2} \]

[Out]

1/2*(10*a*c-3*b^2)/a^2/(-4*a*c+b^2)/e/f^2/(e*x+d)+1/2*(b^2-2*a*c+b*c*(e*x+d)^2)/a/(-4*a*c+b^2)/e/f^2/(e*x+d)/(
a+b*(e*x+d)^2+c*(e*x+d)^4)-1/4*arctan((e*x+d)*2^(1/2)*c^(1/2)/(b-(-4*a*c+b^2)^(1/2))^(1/2))*c^(1/2)*(3*b^3-16*
a*b*c+(-10*a*c+3*b^2)*(-4*a*c+b^2)^(1/2))/a^2/(-4*a*c+b^2)^(3/2)/e/f^2*2^(1/2)/(b-(-4*a*c+b^2)^(1/2))^(1/2)+1/
4*arctan((e*x+d)*2^(1/2)*c^(1/2)/(b+(-4*a*c+b^2)^(1/2))^(1/2))*c^(1/2)*(3*b^3-16*a*b*c-(-10*a*c+3*b^2)*(-4*a*c
+b^2)^(1/2))/a^2/(-4*a*c+b^2)^(3/2)/e/f^2*2^(1/2)/(b+(-4*a*c+b^2)^(1/2))^(1/2)

Rubi [A] (verified)

Time = 1.03 (sec) , antiderivative size = 360, normalized size of antiderivative = 1.00, number of steps used = 6, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.152, Rules used = {1156, 1135, 1295, 1180, 211} \[ \int \frac {1}{(d f+e f x)^2 \left (a+b (d+e x)^2+c (d+e x)^4\right )^2} \, dx=-\frac {\sqrt {c} \left (\left (3 b^2-10 a c\right ) \sqrt {b^2-4 a c}-16 a b c+3 b^3\right ) \arctan \left (\frac {\sqrt {2} \sqrt {c} (d+e x)}{\sqrt {b-\sqrt {b^2-4 a c}}}\right )}{2 \sqrt {2} a^2 e f^2 \left (b^2-4 a c\right )^{3/2} \sqrt {b-\sqrt {b^2-4 a c}}}+\frac {\sqrt {c} \left (-\left (3 b^2-10 a c\right ) \sqrt {b^2-4 a c}-16 a b c+3 b^3\right ) \arctan \left (\frac {\sqrt {2} \sqrt {c} (d+e x)}{\sqrt {\sqrt {b^2-4 a c}+b}}\right )}{2 \sqrt {2} a^2 e f^2 \left (b^2-4 a c\right )^{3/2} \sqrt {\sqrt {b^2-4 a c}+b}}-\frac {3 b^2-10 a c}{2 a^2 e f^2 \left (b^2-4 a c\right ) (d+e x)}+\frac {-2 a c+b^2+b c (d+e x)^2}{2 a e f^2 \left (b^2-4 a c\right ) (d+e x) \left (a+b (d+e x)^2+c (d+e x)^4\right )} \]

[In]

Int[1/((d*f + e*f*x)^2*(a + b*(d + e*x)^2 + c*(d + e*x)^4)^2),x]

[Out]

-1/2*(3*b^2 - 10*a*c)/(a^2*(b^2 - 4*a*c)*e*f^2*(d + e*x)) + (b^2 - 2*a*c + b*c*(d + e*x)^2)/(2*a*(b^2 - 4*a*c)
*e*f^2*(d + e*x)*(a + b*(d + e*x)^2 + c*(d + e*x)^4)) - (Sqrt[c]*(3*b^3 - 16*a*b*c + (3*b^2 - 10*a*c)*Sqrt[b^2
- 4*a*c])*ArcTan[(Sqrt[2]*Sqrt[c]*(d + e*x))/Sqrt[b - Sqrt[b^2 - 4*a*c]]])/(2*Sqrt[2]*a^2*(b^2 - 4*a*c)^(3/2)
*Sqrt[b - Sqrt[b^2 - 4*a*c]]*e*f^2) + (Sqrt[c]*(3*b^3 - 16*a*b*c - (3*b^2 - 10*a*c)*Sqrt[b^2 - 4*a*c])*ArcTan[
(Sqrt[2]*Sqrt[c]*(d + e*x))/Sqrt[b + Sqrt[b^2 - 4*a*c]]])/(2*Sqrt[2]*a^2*(b^2 - 4*a*c)^(3/2)*Sqrt[b + Sqrt[b^2
- 4*a*c]]*e*f^2)

Rule 211

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[a/b, 2]/a)*ArcTan[x/Rt[a/b, 2]], x] /; FreeQ[{a, b}, x]
&& PosQ[a/b]

Rule 1135

Int[((d_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4)^(p_), x_Symbol] :> Simp[(-(d*x)^(m + 1))*(b^2 - 2*
a*c + b*c*x^2)*((a + b*x^2 + c*x^4)^(p + 1)/(2*a*d*(p + 1)*(b^2 - 4*a*c))), x] + Dist[1/(2*a*(p + 1)*(b^2 - 4*
a*c)), Int[(d*x)^m*(a + b*x^2 + c*x^4)^(p + 1)*Simp[b^2*(m + 2*p + 3) - 2*a*c*(m + 4*p + 5) + b*c*(m + 4*p + 7
)*x^2, x], x], x] /; FreeQ[{a, b, c, d, m}, x] && NeQ[b^2 - 4*a*c, 0] && LtQ[p, -1] && IntegerQ[2*p] && (Integ
erQ[p] || IntegerQ[m])

Rule 1156

Int[(u_)^(m_.)*((a_.) + (b_.)*(v_)^2 + (c_.)*(v_)^4)^(p_.), x_Symbol] :> Dist[u^m/(Coefficient[v, x, 1]*v^m),
Subst[Int[x^m*(a + b*x^2 + c*x^(2*2))^p, x], x, v], x] /; FreeQ[{a, b, c, m, p}, x] && LinearPairQ[u, v, x]

Rule 1180

Int[((d_) + (e_.)*(x_)^2)/((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4), x_Symbol] :> With[{q = Rt[b^2 - 4*a*c, 2]}, Di
st[e/2 + (2*c*d - b*e)/(2*q), Int[1/(b/2 - q/2 + c*x^2), x], x] + Dist[e/2 - (2*c*d - b*e)/(2*q), Int[1/(b/2 +
q/2 + c*x^2), x], x]] /; FreeQ[{a, b, c, d, e}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[c*d^2 - a*e^2, 0] && PosQ[b^
2 - 4*a*c]

Rule 1295

Int[((f_.)*(x_))^(m_.)*((d_) + (e_.)*(x_)^2)*((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4)^(p_), x_Symbol] :> Simp[d*(f
*x)^(m + 1)*((a + b*x^2 + c*x^4)^(p + 1)/(a*f*(m + 1))), x] + Dist[1/(a*f^2*(m + 1)), Int[(f*x)^(m + 2)*(a + b
*x^2 + c*x^4)^p*Simp[a*e*(m + 1) - b*d*(m + 2*p + 3) - c*d*(m + 4*p + 5)*x^2, x], x], x] /; FreeQ[{a, b, c, d,
e, f, p}, x] && NeQ[b^2 - 4*a*c, 0] && LtQ[m, -1] && IntegerQ[2*p] && (IntegerQ[p] || IntegerQ[m])

Rubi steps \begin{align*} \text {integral}& = \frac {\text {Subst}\left (\int \frac {1}{x^2 \left (a+b x^2+c x^4\right )^2} \, dx,x,d+e x\right )}{e f^2} \\ & = \frac {b^2-2 a c+b c (d+e x)^2}{2 a \left (b^2-4 a c\right ) e f^2 (d+e x) \left (a+b (d+e x)^2+c (d+e x)^4\right )}-\frac {\text {Subst}\left (\int \frac {-3 b^2+10 a c-3 b c x^2}{x^2 \left (a+b x^2+c x^4\right )} \, dx,x,d+e x\right )}{2 a \left (b^2-4 a c\right ) e f^2} \\ & = -\frac {3 b^2-10 a c}{2 a^2 \left (b^2-4 a c\right ) e f^2 (d+e x)}+\frac {b^2-2 a c+b c (d+e x)^2}{2 a \left (b^2-4 a c\right ) e f^2 (d+e x) \left (a+b (d+e x)^2+c (d+e x)^4\right )}+\frac {\text {Subst}\left (\int \frac {-b \left (3 b^2-13 a c\right )-c \left (3 b^2-10 a c\right ) x^2}{a+b x^2+c x^4} \, dx,x,d+e x\right )}{2 a^2 \left (b^2-4 a c\right ) e f^2} \\ & = -\frac {3 b^2-10 a c}{2 a^2 \left (b^2-4 a c\right ) e f^2 (d+e x)}+\frac {b^2-2 a c+b c (d+e x)^2}{2 a \left (b^2-4 a c\right ) e f^2 (d+e x) \left (a+b (d+e x)^2+c (d+e x)^4\right )}-\frac {\left (c \left (3 b^2-10 a c+\frac {3 b^3}{\sqrt {b^2-4 a c}}-\frac {16 a b c}{\sqrt {b^2-4 a c}}\right )\right ) \text {Subst}\left (\int \frac {1}{\frac {b}{2}-\frac {1}{2} \sqrt {b^2-4 a c}+c x^2} \, dx,x,d+e x\right )}{4 a^2 \left (b^2-4 a c\right ) e f^2}-\frac {\left (c \left (3 b^2-10 a c-\frac {3 b^3}{\sqrt {b^2-4 a c}}+\frac {16 a b c}{\sqrt {b^2-4 a c}}\right )\right ) \text {Subst}\left (\int \frac {1}{\frac {b}{2}+\frac {1}{2} \sqrt {b^2-4 a c}+c x^2} \, dx,x,d+e x\right )}{4 a^2 \left (b^2-4 a c\right ) e f^2} \\ & = -\frac {3 b^2-10 a c}{2 a^2 \left (b^2-4 a c\right ) e f^2 (d+e x)}+\frac {b^2-2 a c+b c (d+e x)^2}{2 a \left (b^2-4 a c\right ) e f^2 (d+e x) \left (a+b (d+e x)^2+c (d+e x)^4\right )}-\frac {\sqrt {c} \left (3 b^2-10 a c+\frac {3 b^3}{\sqrt {b^2-4 a c}}-\frac {16 a b c}{\sqrt {b^2-4 a c}}\right ) \tan ^{-1}\left (\frac {\sqrt {2} \sqrt {c} (d+e x)}{\sqrt {b-\sqrt {b^2-4 a c}}}\right )}{2 \sqrt {2} a^2 \left (b^2-4 a c\right ) \sqrt {b-\sqrt {b^2-4 a c}} e f^2}-\frac {\sqrt {c} \left (3 b^2-10 a c-\frac {3 b^3}{\sqrt {b^2-4 a c}}+\frac {16 a b c}{\sqrt {b^2-4 a c}}\right ) \tan ^{-1}\left (\frac {\sqrt {2} \sqrt {c} (d+e x)}{\sqrt {b+\sqrt {b^2-4 a c}}}\right )}{2 \sqrt {2} a^2 \left (b^2-4 a c\right ) \sqrt {b+\sqrt {b^2-4 a c}} e f^2} \\ \end{align*}

Mathematica [A] (verified)

Time = 0.97 (sec) , antiderivative size = 342, normalized size of antiderivative = 0.95 \[ \int \frac {1}{(d f+e f x)^2 \left (a+b (d+e x)^2+c (d+e x)^4\right )^2} \, dx=\frac {-\frac {4}{d+e x}+\frac {2 (d+e x) \left (b^3-3 a b c+b^2 c (d+e x)^2-2 a c^2 (d+e x)^2\right )}{\left (-b^2+4 a c\right ) \left (a+b (d+e x)^2+c (d+e x)^4\right )}+\frac {\sqrt {2} \sqrt {c} \left (-3 b^3+16 a b c-3 b^2 \sqrt {b^2-4 a c}+10 a c \sqrt {b^2-4 a c}\right ) \arctan \left (\frac {\sqrt {2} \sqrt {c} (d+e x)}{\sqrt {b-\sqrt {b^2-4 a c}}}\right )}{\left (b^2-4 a c\right )^{3/2} \sqrt {b-\sqrt {b^2-4 a c}}}+\frac {\sqrt {2} \sqrt {c} \left (3 b^3-16 a b c-3 b^2 \sqrt {b^2-4 a c}+10 a c \sqrt {b^2-4 a c}\right ) \arctan \left (\frac {\sqrt {2} \sqrt {c} (d+e x)}{\sqrt {b+\sqrt {b^2-4 a c}}}\right )}{\left (b^2-4 a c\right )^{3/2} \sqrt {b+\sqrt {b^2-4 a c}}}}{4 a^2 e f^2} \]

[In]

Integrate[1/((d*f + e*f*x)^2*(a + b*(d + e*x)^2 + c*(d + e*x)^4)^2),x]

[Out]

(-4/(d + e*x) + (2*(d + e*x)*(b^3 - 3*a*b*c + b^2*c*(d + e*x)^2 - 2*a*c^2*(d + e*x)^2))/((-b^2 + 4*a*c)*(a + b
*(d + e*x)^2 + c*(d + e*x)^4)) + (Sqrt[2]*Sqrt[c]*(-3*b^3 + 16*a*b*c - 3*b^2*Sqrt[b^2 - 4*a*c] + 10*a*c*Sqrt[b
^2 - 4*a*c])*ArcTan[(Sqrt[2]*Sqrt[c]*(d + e*x))/Sqrt[b - Sqrt[b^2 - 4*a*c]]])/((b^2 - 4*a*c)^(3/2)*Sqrt[b - Sq
rt[b^2 - 4*a*c]]) + (Sqrt[2]*Sqrt[c]*(3*b^3 - 16*a*b*c - 3*b^2*Sqrt[b^2 - 4*a*c] + 10*a*c*Sqrt[b^2 - 4*a*c])*A
rcTan[(Sqrt[2]*Sqrt[c]*(d + e*x))/Sqrt[b + Sqrt[b^2 - 4*a*c]]])/((b^2 - 4*a*c)^(3/2)*Sqrt[b + Sqrt[b^2 - 4*a*c
]]))/(4*a^2*e*f^2)

Maple [C] (verified)

Result contains higher order function than in optimal. Order 9 vs. order 3.

Time = 0.75 (sec) , antiderivative size = 445, normalized size of antiderivative = 1.24

method	result	size
default	\(\frac {-\frac {\frac {\frac {c \,e^{2} \left (2 a c -b^{2}\right ) x^{3}}{8 a c -2 b^{2}}+\frac {3 d c e \left (2 a c -b^{2}\right ) x^{2}}{2 \left (4 a c -b^{2}\right )}+\frac {\left (6 a \,c^{2} d^{2}-3 b^{2} c \,d^{2}+3 a b c -b^{3}\right ) x}{8 a c -2 b^{2}}+\frac {d \left (2 a \,c^{2} d^{2}-b^{2} c \,d^{2}+3 a b c -b^{3}\right )}{2 e \left (4 a c -b^{2}\right )}}{c \,x^{4} e^{4}+4 c d \,e^{3} x^{3}+6 c \,d^{2} e^{2} x^{2}+4 c \,d^{3} e x +b \,e^{2} x^{2}+d^{4} c +2 b d e x +b \,d^{2}+a}+\frac {\munderset {\textit {\_R} =\operatorname {RootOf}\left (c \,e^{4} \textit {\_Z}^{4}+4 c d \,e^{3} \textit {\_Z}^{3}+\left (6 c \,d^{2} e^{2}+b \,e^{2}\right ) \textit {\_Z}^{2}+\left (4 d^{3} e c +2 b d e \right ) \textit {\_Z} +d^{4} c +b \,d^{2}+a \right )}{\sum }\frac {\left (c \,e^{2} \left (10 a c -3 b^{2}\right ) \textit {\_R}^{2}+2 d c e \left (10 a c -3 b^{2}\right ) \textit {\_R} +10 a \,c^{2} d^{2}-3 b^{2} c \,d^{2}+13 a b c -3 b^{3}\right ) \ln \left (x -\textit {\_R} \right )}{2 e^{3} c \,\textit {\_R}^{3}+6 c d \,e^{2} \textit {\_R}^{2}+6 c \,d^{2} e \textit {\_R} +2 d^{3} c +b e \textit {\_R} +b d}}{4 \left (4 a c -b^{2}\right ) e}}{a^{2}}-\frac {1}{a^{2} e \left (e x +d \right )}}{f^{2}}\)	\(445\)
risch	\(\text {Expression too large to display}\)	\(1268\)

[In]

int(1/(e*f*x+d*f)^2/(a+b*(e*x+d)^2+c*(e*x+d)^4)^2,x,method=_RETURNVERBOSE)

[Out]

1/f^2*(-1/a^2*((1/2*c*e^2*(2*a*c-b^2)/(4*a*c-b^2)*x^3+3/2*d*c*e*(2*a*c-b^2)/(4*a*c-b^2)*x^2+1/2*(6*a*c^2*d^2-3
*b^2*c*d^2+3*a*b*c-b^3)/(4*a*c-b^2)*x+1/2*d/e*(2*a*c^2*d^2-b^2*c*d^2+3*a*b*c-b^3)/(4*a*c-b^2))/(c*e^4*x^4+4*c*
d*e^3*x^3+6*c*d^2*e^2*x^2+4*c*d^3*e*x+b*e^2*x^2+c*d^4+2*b*d*e*x+b*d^2+a)+1/4/(4*a*c-b^2)/e*sum((c*e^2*(10*a*c-
3*b^2)*_R^2+2*d*c*e*(10*a*c-3*b^2)*_R+10*a*c^2*d^2-3*b^2*c*d^2+13*a*b*c-3*b^3)/(2*_R^3*c*e^3+6*_R^2*c*d*e^2+6*
_R*c*d^2*e+2*c*d^3+_R*b*e+b*d)*ln(x-_R),_R=RootOf(c*e^4*_Z^4+4*c*d*e^3*_Z^3+(6*c*d^2*e^2+b*e^2)*_Z^2+(4*c*d^3*
e+2*b*d*e)*_Z+d^4*c+b*d^2+a)))-1/a^2/e/(e*x+d))

Fricas [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 4520 vs. \(2 (312) = 624\).

Time = 0.39 (sec) , antiderivative size = 4520, normalized size of antiderivative = 12.56 \[ \int \frac {1}{(d f+e f x)^2 \left (a+b (d+e x)^2+c (d+e x)^4\right )^2} \, dx=\text {Too large to display} \]

[In]

integrate(1/(e*f*x+d*f)^2/(a+b*(e*x+d)^2+c*(e*x+d)^4)^2,x, algorithm="fricas")

[Out]

-1/4*(2*(3*b^2*c - 10*a*c^2)*e^4*x^4 + 8*(3*b^2*c - 10*a*c^2)*d*e^3*x^3 + 2*(3*b^2*c - 10*a*c^2)*d^4 + 2*(3*b^
3 - 11*a*b*c + 6*(3*b^2*c - 10*a*c^2)*d^2)*e^2*x^2 + 4*a*b^2 - 16*a^2*c + 2*(3*b^3 - 11*a*b*c)*d^2 + 4*(2*(3*b
^2*c - 10*a*c^2)*d^3 + (3*b^3 - 11*a*b*c)*d)*e*x + sqrt(1/2)*((a^2*b^2*c - 4*a^3*c^2)*e^6*f^2*x^5 + 5*(a^2*b^2
*c - 4*a^3*c^2)*d*e^5*f^2*x^4 + (a^2*b^3 - 4*a^3*b*c + 10*(a^2*b^2*c - 4*a^3*c^2)*d^2)*e^4*f^2*x^3 + (10*(a^2*
b^2*c - 4*a^3*c^2)*d^3 + 3*(a^2*b^3 - 4*a^3*b*c)*d)*e^3*f^2*x^2 + (a^3*b^2 - 4*a^4*c + 5*(a^2*b^2*c - 4*a^3*c^
2)*d^4 + 3*(a^2*b^3 - 4*a^3*b*c)*d^2)*e^2*f^2*x + ((a^2*b^2*c - 4*a^3*c^2)*d^5 + (a^2*b^3 - 4*a^3*b*c)*d^3 + (
a^3*b^2 - 4*a^4*c)*d)*e*f^2)*sqrt(-((a^5*b^6 - 12*a^6*b^4*c + 48*a^7*b^2*c^2 - 64*a^8*c^3)*e^2*f^4*sqrt((81*b^
8 - 918*a*b^6*c + 3051*a^2*b^4*c^2 - 2550*a^3*b^2*c^3 + 625*a^4*c^4)/((a^10*b^6 - 12*a^11*b^4*c + 48*a^12*b^2*
c^2 - 64*a^13*c^3)*e^4*f^8)) + 9*b^7 - 105*a*b^5*c + 385*a^2*b^3*c^2 - 420*a^3*b*c^3)/((a^5*b^6 - 12*a^6*b^4*c
+ 48*a^7*b^2*c^2 - 64*a^8*c^3)*e^2*f^4))*log(-(189*b^6*c^3 - 1971*a*b^4*c^4 + 5625*a^2*b^2*c^5 - 2500*a^3*c^6
)*e*x - (189*b^6*c^3 - 1971*a*b^4*c^4 + 5625*a^2*b^2*c^5 - 2500*a^3*c^6)*d + 1/2*sqrt(1/2)*((3*a^5*b^10 - 55*a
^6*b^8*c + 392*a^7*b^6*c^2 - 1344*a^8*b^4*c^3 + 2176*a^9*b^2*c^4 - 1280*a^10*c^5)*e^3*f^6*sqrt((81*b^8 - 918*a
*b^6*c + 3051*a^2*b^4*c^2 - 2550*a^3*b^2*c^3 + 625*a^4*c^4)/((a^10*b^6 - 12*a^11*b^4*c + 48*a^12*b^2*c^2 - 64*
a^13*c^3)*e^4*f^8)) - (27*b^11 - 486*a*b^9*c + 3330*a^2*b^7*c^2 - 10549*a^3*b^5*c^3 + 14408*a^4*b^3*c^4 - 5200
*a^5*b*c^5)*e*f^2)*sqrt(-((a^5*b^6 - 12*a^6*b^4*c + 48*a^7*b^2*c^2 - 64*a^8*c^3)*e^2*f^4*sqrt((81*b^8 - 918*a*
b^6*c + 3051*a^2*b^4*c^2 - 2550*a^3*b^2*c^3 + 625*a^4*c^4)/((a^10*b^6 - 12*a^11*b^4*c + 48*a^12*b^2*c^2 - 64*a
^13*c^3)*e^4*f^8)) + 9*b^7 - 105*a*b^5*c + 385*a^2*b^3*c^2 - 420*a^3*b*c^3)/((a^5*b^6 - 12*a^6*b^4*c + 48*a^7*
b^2*c^2 - 64*a^8*c^3)*e^2*f^4))) - sqrt(1/2)*((a^2*b^2*c - 4*a^3*c^2)*e^6*f^2*x^5 + 5*(a^2*b^2*c - 4*a^3*c^2)*
d*e^5*f^2*x^4 + (a^2*b^3 - 4*a^3*b*c + 10*(a^2*b^2*c - 4*a^3*c^2)*d^2)*e^4*f^2*x^3 + (10*(a^2*b^2*c - 4*a^3*c^
2)*d^3 + 3*(a^2*b^3 - 4*a^3*b*c)*d)*e^3*f^2*x^2 + (a^3*b^2 - 4*a^4*c + 5*(a^2*b^2*c - 4*a^3*c^2)*d^4 + 3*(a^2*
b^3 - 4*a^3*b*c)*d^2)*e^2*f^2*x + ((a^2*b^2*c - 4*a^3*c^2)*d^5 + (a^2*b^3 - 4*a^3*b*c)*d^3 + (a^3*b^2 - 4*a^4*
c)*d)*e*f^2)*sqrt(-((a^5*b^6 - 12*a^6*b^4*c + 48*a^7*b^2*c^2 - 64*a^8*c^3)*e^2*f^4*sqrt((81*b^8 - 918*a*b^6*c
+ 3051*a^2*b^4*c^2 - 2550*a^3*b^2*c^3 + 625*a^4*c^4)/((a^10*b^6 - 12*a^11*b^4*c + 48*a^12*b^2*c^2 - 64*a^13*c^
3)*e^4*f^8)) + 9*b^7 - 105*a*b^5*c + 385*a^2*b^3*c^2 - 420*a^3*b*c^3)/((a^5*b^6 - 12*a^6*b^4*c + 48*a^7*b^2*c^
2 - 64*a^8*c^3)*e^2*f^4))*log(-(189*b^6*c^3 - 1971*a*b^4*c^4 + 5625*a^2*b^2*c^5 - 2500*a^3*c^6)*e*x - (189*b^6
*c^3 - 1971*a*b^4*c^4 + 5625*a^2*b^2*c^5 - 2500*a^3*c^6)*d - 1/2*sqrt(1/2)*((3*a^5*b^10 - 55*a^6*b^8*c + 392*a
^7*b^6*c^2 - 1344*a^8*b^4*c^3 + 2176*a^9*b^2*c^4 - 1280*a^10*c^5)*e^3*f^6*sqrt((81*b^8 - 918*a*b^6*c + 3051*a^
2*b^4*c^2 - 2550*a^3*b^2*c^3 + 625*a^4*c^4)/((a^10*b^6 - 12*a^11*b^4*c + 48*a^12*b^2*c^2 - 64*a^13*c^3)*e^4*f^
8)) - (27*b^11 - 486*a*b^9*c + 3330*a^2*b^7*c^2 - 10549*a^3*b^5*c^3 + 14408*a^4*b^3*c^4 - 5200*a^5*b*c^5)*e*f^
2)*sqrt(-((a^5*b^6 - 12*a^6*b^4*c + 48*a^7*b^2*c^2 - 64*a^8*c^3)*e^2*f^4*sqrt((81*b^8 - 918*a*b^6*c + 3051*a^2
*b^4*c^2 - 2550*a^3*b^2*c^3 + 625*a^4*c^4)/((a^10*b^6 - 12*a^11*b^4*c + 48*a^12*b^2*c^2 - 64*a^13*c^3)*e^4*f^8
)) + 9*b^7 - 105*a*b^5*c + 385*a^2*b^3*c^2 - 420*a^3*b*c^3)/((a^5*b^6 - 12*a^6*b^4*c + 48*a^7*b^2*c^2 - 64*a^8
*c^3)*e^2*f^4))) - sqrt(1/2)*((a^2*b^2*c - 4*a^3*c^2)*e^6*f^2*x^5 + 5*(a^2*b^2*c - 4*a^3*c^2)*d*e^5*f^2*x^4 +
(a^2*b^3 - 4*a^3*b*c + 10*(a^2*b^2*c - 4*a^3*c^2)*d^2)*e^4*f^2*x^3 + (10*(a^2*b^2*c - 4*a^3*c^2)*d^3 + 3*(a^2*
b^3 - 4*a^3*b*c)*d)*e^3*f^2*x^2 + (a^3*b^2 - 4*a^4*c + 5*(a^2*b^2*c - 4*a^3*c^2)*d^4 + 3*(a^2*b^3 - 4*a^3*b*c)
*d^2)*e^2*f^2*x + ((a^2*b^2*c - 4*a^3*c^2)*d^5 + (a^2*b^3 - 4*a^3*b*c)*d^3 + (a^3*b^2 - 4*a^4*c)*d)*e*f^2)*sqr
t(((a^5*b^6 - 12*a^6*b^4*c + 48*a^7*b^2*c^2 - 64*a^8*c^3)*e^2*f^4*sqrt((81*b^8 - 918*a*b^6*c + 3051*a^2*b^4*c^
2 - 2550*a^3*b^2*c^3 + 625*a^4*c^4)/((a^10*b^6 - 12*a^11*b^4*c + 48*a^12*b^2*c^2 - 64*a^13*c^3)*e^4*f^8)) - 9*
b^7 + 105*a*b^5*c - 385*a^2*b^3*c^2 + 420*a^3*b*c^3)/((a^5*b^6 - 12*a^6*b^4*c + 48*a^7*b^2*c^2 - 64*a^8*c^3)*e
^2*f^4))*log(-(189*b^6*c^3 - 1971*a*b^4*c^4 + 5625*a^2*b^2*c^5 - 2500*a^3*c^6)*e*x - (189*b^6*c^3 - 1971*a*b^4
*c^4 + 5625*a^2*b^2*c^5 - 2500*a^3*c^6)*d + 1/2*sqrt(1/2)*((3*a^5*b^10 - 55*a^6*b^8*c + 392*a^7*b^6*c^2 - 1344
*a^8*b^4*c^3 + 2176*a^9*b^2*c^4 - 1280*a^10*c^5)*e^3*f^6*sqrt((81*b^8 - 918*a*b^6*c + 3051*a^2*b^4*c^2 - 2550*
a^3*b^2*c^3 + 625*a^4*c^4)/((a^10*b^6 - 12*a^11*b^4*c + 48*a^12*b^2*c^2 - 64*a^13*c^3)*e^4*f^8)) + (27*b^11 -
486*a*b^9*c + 3330*a^2*b^7*c^2 - 10549*a^3*b^5*c^3 + 14408*a^4*b^3*c^4 - 5200*a^5*b*c^5)*e*f^2)*sqrt(((a^5*b^6
- 12*a^6*b^4*c + 48*a^7*b^2*c^2 - 64*a^8*c^3)*e^2*f^4*sqrt((81*b^8 - 918*a*b^6*c + 3051*a^2*b^4*c^2 - 2550*a^
3*b^2*c^3 + 625*a^4*c^4)/((a^10*b^6 - 12*a^11*b^4*c + 48*a^12*b^2*c^2 - 64*a^13*c^3)*e^4*f^8)) - 9*b^7 + 105*a
*b^5*c - 385*a^2*b^3*c^2 + 420*a^3*b*c^3)/((a^5*b^6 - 12*a^6*b^4*c + 48*a^7*b^2*c^2 - 64*a^8*c^3)*e^2*f^4))) +
sqrt(1/2)*((a^2*b^2*c - 4*a^3*c^2)*e^6*f^2*x^5 + 5*(a^2*b^2*c - 4*a^3*c^2)*d*e^5*f^2*x^4 + (a^2*b^3 - 4*a^3*b
*c + 10*(a^2*b^2*c - 4*a^3*c^2)*d^2)*e^4*f^2*x^3 + (10*(a^2*b^2*c - 4*a^3*c^2)*d^3 + 3*(a^2*b^3 - 4*a^3*b*c)*d
)*e^3*f^2*x^2 + (a^3*b^2 - 4*a^4*c + 5*(a^2*b^2*c - 4*a^3*c^2)*d^4 + 3*(a^2*b^3 - 4*a^3*b*c)*d^2)*e^2*f^2*x +
((a^2*b^2*c - 4*a^3*c^2)*d^5 + (a^2*b^3 - 4*a^3*b*c)*d^3 + (a^3*b^2 - 4*a^4*c)*d)*e*f^2)*sqrt(((a^5*b^6 - 12*a
^6*b^4*c + 48*a^7*b^2*c^2 - 64*a^8*c^3)*e^2*f^4*sqrt((81*b^8 - 918*a*b^6*c + 3051*a^2*b^4*c^2 - 2550*a^3*b^2*c
^3 + 625*a^4*c^4)/((a^10*b^6 - 12*a^11*b^4*c + 48*a^12*b^2*c^2 - 64*a^13*c^3)*e^4*f^8)) - 9*b^7 + 105*a*b^5*c
- 385*a^2*b^3*c^2 + 420*a^3*b*c^3)/((a^5*b^6 - 12*a^6*b^4*c + 48*a^7*b^2*c^2 - 64*a^8*c^3)*e^2*f^4))*log(-(189
*b^6*c^3 - 1971*a*b^4*c^4 + 5625*a^2*b^2*c^5 - 2500*a^3*c^6)*e*x - (189*b^6*c^3 - 1971*a*b^4*c^4 + 5625*a^2*b^
2*c^5 - 2500*a^3*c^6)*d - 1/2*sqrt(1/2)*((3*a^5*b^10 - 55*a^6*b^8*c + 392*a^7*b^6*c^2 - 1344*a^8*b^4*c^3 + 217
6*a^9*b^2*c^4 - 1280*a^10*c^5)*e^3*f^6*sqrt((81*b^8 - 918*a*b^6*c + 3051*a^2*b^4*c^2 - 2550*a^3*b^2*c^3 + 625*
a^4*c^4)/((a^10*b^6 - 12*a^11*b^4*c + 48*a^12*b^2*c^2 - 64*a^13*c^3)*e^4*f^8)) + (27*b^11 - 486*a*b^9*c + 3330
*a^2*b^7*c^2 - 10549*a^3*b^5*c^3 + 14408*a^4*b^3*c^4 - 5200*a^5*b*c^5)*e*f^2)*sqrt(((a^5*b^6 - 12*a^6*b^4*c +
48*a^7*b^2*c^2 - 64*a^8*c^3)*e^2*f^4*sqrt((81*b^8 - 918*a*b^6*c + 3051*a^2*b^4*c^2 - 2550*a^3*b^2*c^3 + 625*a^
4*c^4)/((a^10*b^6 - 12*a^11*b^4*c + 48*a^12*b^2*c^2 - 64*a^13*c^3)*e^4*f^8)) - 9*b^7 + 105*a*b^5*c - 385*a^2*b
^3*c^2 + 420*a^3*b*c^3)/((a^5*b^6 - 12*a^6*b^4*c + 48*a^7*b^2*c^2 - 64*a^8*c^3)*e^2*f^4))))/((a^2*b^2*c - 4*a^
3*c^2)*e^6*f^2*x^5 + 5*(a^2*b^2*c - 4*a^3*c^2)*d*e^5*f^2*x^4 + (a^2*b^3 - 4*a^3*b*c + 10*(a^2*b^2*c - 4*a^3*c^
2)*d^2)*e^4*f^2*x^3 + (10*(a^2*b^2*c - 4*a^3*c^2)*d^3 + 3*(a^2*b^3 - 4*a^3*b*c)*d)*e^3*f^2*x^2 + (a^3*b^2 - 4*
a^4*c + 5*(a^2*b^2*c - 4*a^3*c^2)*d^4 + 3*(a^2*b^3 - 4*a^3*b*c)*d^2)*e^2*f^2*x + ((a^2*b^2*c - 4*a^3*c^2)*d^5
+ (a^2*b^3 - 4*a^3*b*c)*d^3 + (a^3*b^2 - 4*a^4*c)*d)*e*f^2)

Sympy [F(-1)]

Timed out. \[ \int \frac {1}{(d f+e f x)^2 \left (a+b (d+e x)^2+c (d+e x)^4\right )^2} \, dx=\text {Timed out} \]

[In]

integrate(1/(e*f*x+d*f)**2/(a+b*(e*x+d)**2+c*(e*x+d)**4)**2,x)

[Out]

Timed out

Maxima [F]

\[ \int \frac {1}{(d f+e f x)^2 \left (a+b (d+e x)^2+c (d+e x)^4\right )^2} \, dx=\int { \frac {1}{{\left ({\left (e x + d\right )}^{4} c + {\left (e x + d\right )}^{2} b + a\right )}^{2} {\left (e f x + d f\right )}^{2}} \,d x } \]

[In]

integrate(1/(e*f*x+d*f)^2/(a+b*(e*x+d)^2+c*(e*x+d)^4)^2,x, algorithm="maxima")

[Out]

-1/2*((3*b^2*c - 10*a*c^2)*e^4*x^4 + 4*(3*b^2*c - 10*a*c^2)*d*e^3*x^3 + (3*b^2*c - 10*a*c^2)*d^4 + (3*b^3 - 11
*a*b*c + 6*(3*b^2*c - 10*a*c^2)*d^2)*e^2*x^2 + 2*a*b^2 - 8*a^2*c + (3*b^3 - 11*a*b*c)*d^2 + 2*(2*(3*b^2*c - 10
*a*c^2)*d^3 + (3*b^3 - 11*a*b*c)*d)*e*x)/((a^2*b^2*c - 4*a^3*c^2)*e^6*f^2*x^5 + 5*(a^2*b^2*c - 4*a^3*c^2)*d*e^
5*f^2*x^4 + (a^2*b^3 - 4*a^3*b*c + 10*(a^2*b^2*c - 4*a^3*c^2)*d^2)*e^4*f^2*x^3 + (10*(a^2*b^2*c - 4*a^3*c^2)*d
^3 + 3*(a^2*b^3 - 4*a^3*b*c)*d)*e^3*f^2*x^2 + (a^3*b^2 - 4*a^4*c + 5*(a^2*b^2*c - 4*a^3*c^2)*d^4 + 3*(a^2*b^3
- 4*a^3*b*c)*d^2)*e^2*f^2*x + ((a^2*b^2*c - 4*a^3*c^2)*d^5 + (a^2*b^3 - 4*a^3*b*c)*d^3 + (a^3*b^2 - 4*a^4*c)*d
)*e*f^2) - 1/2*integrate(((3*b^2*c - 10*a*c^2)*e^2*x^2 + 2*(3*b^2*c - 10*a*c^2)*d*e*x + 3*b^3 - 13*a*b*c + (3*
b^2*c - 10*a*c^2)*d^2)/((b^2*c - 4*a*c^2)*e^4*x^4 + 4*(b^2*c - 4*a*c^2)*d*e^3*x^3 + (b^2*c - 4*a*c^2)*d^4 + (b
^3 - 4*a*b*c + 6*(b^2*c - 4*a*c^2)*d^2)*e^2*x^2 + a*b^2 - 4*a^2*c + (b^3 - 4*a*b*c)*d^2 + 2*(2*(b^2*c - 4*a*c^
2)*d^3 + (b^3 - 4*a*b*c)*d)*e*x), x)/(a^2*f^2)

Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 1031 vs. \(2 (312) = 624\).

Time = 0.34 (sec) , antiderivative size = 1031, normalized size of antiderivative = 2.86 \[ \int \frac {1}{(d f+e f x)^2 \left (a+b (d+e x)^2+c (d+e x)^4\right )^2} \, dx=-\frac {\frac {b^{2} c}{{\left (e f x + d f\right )} e f} - \frac {2 \, a c^{2}}{{\left (e f x + d f\right )} e f} + \frac {b^{3} f}{{\left (e f x + d f\right )}^{3} e} - \frac {3 \, a b c f}{{\left (e f x + d f\right )}^{3} e}}{2 \, {\left (a^{2} b^{2} - 4 \, a^{3} c\right )} {\left (c + \frac {b f^{2}}{{\left (e f x + d f\right )}^{2}} + \frac {a f^{4}}{{\left (e f x + d f\right )}^{4}}\right )}} - \frac {1}{{\left (e f x + d f\right )} a^{2} e f} + \frac {{\left ({\left (3 \, a^{4} b^{7} - 31 \, a^{5} b^{5} c + 96 \, a^{6} b^{3} c^{2} - 80 \, a^{7} b c^{3}\right )} \sqrt {2 \, a b + 2 \, \sqrt {b^{2} - 4 \, a c} a} e^{4} f^{8} + 2 \, {\left (3 \, a^{3} b^{2} c - 10 \, a^{4} c^{2}\right )} \sqrt {2 \, a b + 2 \, \sqrt {b^{2} - 4 \, a c} a} \sqrt {b^{2} - 4 \, a c} e^{2} f^{4} {\left | a^{2} b^{2} e^{2} f^{4} - 4 \, a^{3} c e^{2} f^{4} \right |} - {\left (a^{2} b^{2} e^{2} f^{4} - 4 \, a^{3} c e^{2} f^{4}\right )}^{2} {\left (3 \, b^{3} - 13 \, a b c\right )} \sqrt {2 \, a b + 2 \, \sqrt {b^{2} - 4 \, a c} a}\right )} \arctan \left (\frac {2 \, \sqrt {\frac {1}{2}}}{{\left (e f x + d f\right )} e f \sqrt {\frac {a^{2} b^{3} e^{2} f^{4} - 4 \, a^{3} b c e^{2} f^{4} + \sqrt {{\left (a^{2} b^{3} e^{2} f^{4} - 4 \, a^{3} b c e^{2} f^{4}\right )}^{2} - 4 \, {\left (a^{3} b^{2} e^{4} f^{8} - 4 \, a^{4} c e^{4} f^{8}\right )} {\left (a^{2} b^{2} c - 4 \, a^{3} c^{2}\right )}}}{a^{3} b^{2} e^{4} f^{8} - 4 \, a^{4} c e^{4} f^{8}}}}\right )}{16 \, {\left (a^{5} b^{2} c - 4 \, a^{6} c^{2}\right )} \sqrt {b^{2} - 4 \, a c} e^{3} f^{6} {\left | a^{2} b^{2} e^{2} f^{4} - 4 \, a^{3} c e^{2} f^{4} \right |} {\left | a \right |}} - \frac {{\left ({\left (3 \, a^{4} b^{7} - 31 \, a^{5} b^{5} c + 96 \, a^{6} b^{3} c^{2} - 80 \, a^{7} b c^{3}\right )} \sqrt {2 \, a b - 2 \, \sqrt {b^{2} - 4 \, a c} a} e^{4} f^{8} - 2 \, {\left (3 \, a^{3} b^{2} c - 10 \, a^{4} c^{2}\right )} \sqrt {2 \, a b - 2 \, \sqrt {b^{2} - 4 \, a c} a} \sqrt {b^{2} - 4 \, a c} e^{2} f^{4} {\left | a^{2} b^{2} e^{2} f^{4} - 4 \, a^{3} c e^{2} f^{4} \right |} - {\left (a^{2} b^{2} e^{2} f^{4} - 4 \, a^{3} c e^{2} f^{4}\right )}^{2} {\left (3 \, b^{3} - 13 \, a b c\right )} \sqrt {2 \, a b - 2 \, \sqrt {b^{2} - 4 \, a c} a}\right )} \arctan \left (\frac {2 \, \sqrt {\frac {1}{2}}}{{\left (e f x + d f\right )} e f \sqrt {\frac {a^{2} b^{3} e^{2} f^{4} - 4 \, a^{3} b c e^{2} f^{4} - \sqrt {{\left (a^{2} b^{3} e^{2} f^{4} - 4 \, a^{3} b c e^{2} f^{4}\right )}^{2} - 4 \, {\left (a^{3} b^{2} e^{4} f^{8} - 4 \, a^{4} c e^{4} f^{8}\right )} {\left (a^{2} b^{2} c - 4 \, a^{3} c^{2}\right )}}}{a^{3} b^{2} e^{4} f^{8} - 4 \, a^{4} c e^{4} f^{8}}}}\right )}{16 \, {\left (a^{5} b^{2} c - 4 \, a^{6} c^{2}\right )} \sqrt {b^{2} - 4 \, a c} e^{3} f^{6} {\left | a^{2} b^{2} e^{2} f^{4} - 4 \, a^{3} c e^{2} f^{4} \right |} {\left | a \right |}} \]

[In]

integrate(1/(e*f*x+d*f)^2/(a+b*(e*x+d)^2+c*(e*x+d)^4)^2,x, algorithm="giac")

[Out]

-1/2*(b^2*c/((e*f*x + d*f)*e*f) - 2*a*c^2/((e*f*x + d*f)*e*f) + b^3*f/((e*f*x + d*f)^3*e) - 3*a*b*c*f/((e*f*x
+ d*f)^3*e))/((a^2*b^2 - 4*a^3*c)*(c + b*f^2/(e*f*x + d*f)^2 + a*f^4/(e*f*x + d*f)^4)) - 1/((e*f*x + d*f)*a^2*
e*f) + 1/16*((3*a^4*b^7 - 31*a^5*b^5*c + 96*a^6*b^3*c^2 - 80*a^7*b*c^3)*sqrt(2*a*b + 2*sqrt(b^2 - 4*a*c)*a)*e^
4*f^8 + 2*(3*a^3*b^2*c - 10*a^4*c^2)*sqrt(2*a*b + 2*sqrt(b^2 - 4*a*c)*a)*sqrt(b^2 - 4*a*c)*e^2*f^4*abs(a^2*b^2
*e^2*f^4 - 4*a^3*c*e^2*f^4) - (a^2*b^2*e^2*f^4 - 4*a^3*c*e^2*f^4)^2*(3*b^3 - 13*a*b*c)*sqrt(2*a*b + 2*sqrt(b^2
- 4*a*c)*a))*arctan(2*sqrt(1/2)/((e*f*x + d*f)*e*f*sqrt((a^2*b^3*e^2*f^4 - 4*a^3*b*c*e^2*f^4 + sqrt((a^2*b^3*
e^2*f^4 - 4*a^3*b*c*e^2*f^4)^2 - 4*(a^3*b^2*e^4*f^8 - 4*a^4*c*e^4*f^8)*(a^2*b^2*c - 4*a^3*c^2)))/(a^3*b^2*e^4*
f^8 - 4*a^4*c*e^4*f^8))))/((a^5*b^2*c - 4*a^6*c^2)*sqrt(b^2 - 4*a*c)*e^3*f^6*abs(a^2*b^2*e^2*f^4 - 4*a^3*c*e^2
*f^4)*abs(a)) - 1/16*((3*a^4*b^7 - 31*a^5*b^5*c + 96*a^6*b^3*c^2 - 80*a^7*b*c^3)*sqrt(2*a*b - 2*sqrt(b^2 - 4*a
*c)*a)*e^4*f^8 - 2*(3*a^3*b^2*c - 10*a^4*c^2)*sqrt(2*a*b - 2*sqrt(b^2 - 4*a*c)*a)*sqrt(b^2 - 4*a*c)*e^2*f^4*ab
s(a^2*b^2*e^2*f^4 - 4*a^3*c*e^2*f^4) - (a^2*b^2*e^2*f^4 - 4*a^3*c*e^2*f^4)^2*(3*b^3 - 13*a*b*c)*sqrt(2*a*b - 2
*sqrt(b^2 - 4*a*c)*a))*arctan(2*sqrt(1/2)/((e*f*x + d*f)*e*f*sqrt((a^2*b^3*e^2*f^4 - 4*a^3*b*c*e^2*f^4 - sqrt(
(a^2*b^3*e^2*f^4 - 4*a^3*b*c*e^2*f^4)^2 - 4*(a^3*b^2*e^4*f^8 - 4*a^4*c*e^4*f^8)*(a^2*b^2*c - 4*a^3*c^2)))/(a^3
*b^2*e^4*f^8 - 4*a^4*c*e^4*f^8))))/((a^5*b^2*c - 4*a^6*c^2)*sqrt(b^2 - 4*a*c)*e^3*f^6*abs(a^2*b^2*e^2*f^4 - 4*
a^3*c*e^2*f^4)*abs(a))

Mupad [B] (verification not implemented)

Time = 11.87 (sec) , antiderivative size = 12008, normalized size of antiderivative = 33.36 \[ \int \frac {1}{(d f+e f x)^2 \left (a+b (d+e x)^2+c (d+e x)^4\right )^2} \, dx=\text {Too large to display} \]

[In]

int(1/((d*f + e*f*x)^2*(a + b*(d + e*x)^2 + c*(d + e*x)^4)^2),x)

[Out]

- atan(((-(9*b^13 - 9*b^4*(-(4*a*c - b^2)^9)^(1/2) + 26880*a^6*b*c^6 + 2077*a^2*b^9*c^2 - 10656*a^3*b^7*c^3 +
30240*a^4*b^5*c^4 - 44800*a^5*b^3*c^5 - 25*a^2*c^2*(-(4*a*c - b^2)^9)^(1/2) - 213*a*b^11*c + 51*a*b^2*c*(-(4*a
*c - b^2)^9)^(1/2))/(32*(a^5*b^12*e^2*f^4 + 4096*a^11*c^6*e^2*f^4 + 240*a^7*b^8*c^2*e^2*f^4 - 1280*a^8*b^6*c^3
*e^2*f^4 + 3840*a^9*b^4*c^4*e^2*f^4 - 6144*a^10*b^2*c^5*e^2*f^4 - 24*a^6*b^10*c*e^2*f^4)))^(1/2)*((-(9*b^13 -
9*b^4*(-(4*a*c - b^2)^9)^(1/2) + 26880*a^6*b*c^6 + 2077*a^2*b^9*c^2 - 10656*a^3*b^7*c^3 + 30240*a^4*b^5*c^4 -
44800*a^5*b^3*c^5 - 25*a^2*c^2*(-(4*a*c - b^2)^9)^(1/2) - 213*a*b^11*c + 51*a*b^2*c*(-(4*a*c - b^2)^9)^(1/2))/
(32*(a^5*b^12*e^2*f^4 + 4096*a^11*c^6*e^2*f^4 + 240*a^7*b^8*c^2*e^2*f^4 - 1280*a^8*b^6*c^3*e^2*f^4 + 3840*a^9*
b^4*c^4*e^2*f^4 - 6144*a^10*b^2*c^5*e^2*f^4 - 24*a^6*b^10*c*e^2*f^4)))^(1/2)*((-(9*b^13 - 9*b^4*(-(4*a*c - b^2
)^9)^(1/2) + 26880*a^6*b*c^6 + 2077*a^2*b^9*c^2 - 10656*a^3*b^7*c^3 + 30240*a^4*b^5*c^4 - 44800*a^5*b^3*c^5 -
25*a^2*c^2*(-(4*a*c - b^2)^9)^(1/2) - 213*a*b^11*c + 51*a*b^2*c*(-(4*a*c - b^2)^9)^(1/2))/(32*(a^5*b^12*e^2*f^
4 + 4096*a^11*c^6*e^2*f^4 + 240*a^7*b^8*c^2*e^2*f^4 - 1280*a^8*b^6*c^3*e^2*f^4 + 3840*a^9*b^4*c^4*e^2*f^4 - 61
44*a^10*b^2*c^5*e^2*f^4 - 24*a^6*b^10*c*e^2*f^4)))^(1/2)*(x*(256*a^10*b^13*c^2*e^14*f^10 - 6144*a^11*b^11*c^3*
e^14*f^10 + 61440*a^12*b^9*c^4*e^14*f^10 - 327680*a^13*b^7*c^5*e^14*f^10 + 983040*a^14*b^5*c^6*e^14*f^10 - 157
2864*a^15*b^3*c^7*e^14*f^10 + 1048576*a^16*b*c^8*e^14*f^10) + 1048576*a^16*b*c^8*d*e^13*f^10 + 256*a^10*b^13*c
^2*d*e^13*f^10 - 6144*a^11*b^11*c^3*d*e^13*f^10 + 61440*a^12*b^9*c^4*d*e^13*f^10 - 327680*a^13*b^7*c^5*d*e^13*
f^10 + 983040*a^14*b^5*c^6*d*e^13*f^10 - 1572864*a^15*b^3*c^7*d*e^13*f^10) - 192*a^8*b^13*c^2*e^12*f^8 + 4672*
a^9*b^11*c^3*e^12*f^8 - 47360*a^10*b^9*c^4*e^12*f^8 + 256000*a^11*b^7*c^5*e^12*f^8 - 778240*a^12*b^5*c^6*e^12*
f^8 + 1261568*a^13*b^3*c^7*e^12*f^8 - 851968*a^14*b*c^8*e^12*f^8) + x*(204800*a^12*c^9*e^12*f^6 + 144*a^6*b^12
*c^3*e^12*f^6 - 3264*a^7*b^10*c^4*e^12*f^6 + 30112*a^8*b^8*c^5*e^12*f^6 - 143360*a^9*b^6*c^6*e^12*f^6 + 365568
*a^10*b^4*c^7*e^12*f^6 - 458752*a^11*b^2*c^8*e^12*f^6) + 204800*a^12*c^9*d*e^11*f^6 + 144*a^6*b^12*c^3*d*e^11*
f^6 - 3264*a^7*b^10*c^4*d*e^11*f^6 + 30112*a^8*b^8*c^5*d*e^11*f^6 - 143360*a^9*b^6*c^6*d*e^11*f^6 + 365568*a^1
0*b^4*c^7*d*e^11*f^6 - 458752*a^11*b^2*c^8*d*e^11*f^6)*1i + (-(9*b^13 - 9*b^4*(-(4*a*c - b^2)^9)^(1/2) + 26880
*a^6*b*c^6 + 2077*a^2*b^9*c^2 - 10656*a^3*b^7*c^3 + 30240*a^4*b^5*c^4 - 44800*a^5*b^3*c^5 - 25*a^2*c^2*(-(4*a*
c - b^2)^9)^(1/2) - 213*a*b^11*c + 51*a*b^2*c*(-(4*a*c - b^2)^9)^(1/2))/(32*(a^5*b^12*e^2*f^4 + 4096*a^11*c^6*
e^2*f^4 + 240*a^7*b^8*c^2*e^2*f^4 - 1280*a^8*b^6*c^3*e^2*f^4 + 3840*a^9*b^4*c^4*e^2*f^4 - 6144*a^10*b^2*c^5*e^
2*f^4 - 24*a^6*b^10*c*e^2*f^4)))^(1/2)*((-(9*b^13 - 9*b^4*(-(4*a*c - b^2)^9)^(1/2) + 26880*a^6*b*c^6 + 2077*a^
2*b^9*c^2 - 10656*a^3*b^7*c^3 + 30240*a^4*b^5*c^4 - 44800*a^5*b^3*c^5 - 25*a^2*c^2*(-(4*a*c - b^2)^9)^(1/2) -
213*a*b^11*c + 51*a*b^2*c*(-(4*a*c - b^2)^9)^(1/2))/(32*(a^5*b^12*e^2*f^4 + 4096*a^11*c^6*e^2*f^4 + 240*a^7*b^
8*c^2*e^2*f^4 - 1280*a^8*b^6*c^3*e^2*f^4 + 3840*a^9*b^4*c^4*e^2*f^4 - 6144*a^10*b^2*c^5*e^2*f^4 - 24*a^6*b^10*
c*e^2*f^4)))^(1/2)*((-(9*b^13 - 9*b^4*(-(4*a*c - b^2)^9)^(1/2) + 26880*a^6*b*c^6 + 2077*a^2*b^9*c^2 - 10656*a^
3*b^7*c^3 + 30240*a^4*b^5*c^4 - 44800*a^5*b^3*c^5 - 25*a^2*c^2*(-(4*a*c - b^2)^9)^(1/2) - 213*a*b^11*c + 51*a*
b^2*c*(-(4*a*c - b^2)^9)^(1/2))/(32*(a^5*b^12*e^2*f^4 + 4096*a^11*c^6*e^2*f^4 + 240*a^7*b^8*c^2*e^2*f^4 - 1280
*a^8*b^6*c^3*e^2*f^4 + 3840*a^9*b^4*c^4*e^2*f^4 - 6144*a^10*b^2*c^5*e^2*f^4 - 24*a^6*b^10*c*e^2*f^4)))^(1/2)*(
x*(256*a^10*b^13*c^2*e^14*f^10 - 6144*a^11*b^11*c^3*e^14*f^10 + 61440*a^12*b^9*c^4*e^14*f^10 - 327680*a^13*b^7
*c^5*e^14*f^10 + 983040*a^14*b^5*c^6*e^14*f^10 - 1572864*a^15*b^3*c^7*e^14*f^10 + 1048576*a^16*b*c^8*e^14*f^10
) + 1048576*a^16*b*c^8*d*e^13*f^10 + 256*a^10*b^13*c^2*d*e^13*f^10 - 6144*a^11*b^11*c^3*d*e^13*f^10 + 61440*a^
12*b^9*c^4*d*e^13*f^10 - 327680*a^13*b^7*c^5*d*e^13*f^10 + 983040*a^14*b^5*c^6*d*e^13*f^10 - 1572864*a^15*b^3*
c^7*d*e^13*f^10) + 192*a^8*b^13*c^2*e^12*f^8 - 4672*a^9*b^11*c^3*e^12*f^8 + 47360*a^10*b^9*c^4*e^12*f^8 - 2560
00*a^11*b^7*c^5*e^12*f^8 + 778240*a^12*b^5*c^6*e^12*f^8 - 1261568*a^13*b^3*c^7*e^12*f^8 + 851968*a^14*b*c^8*e^
12*f^8) + x*(204800*a^12*c^9*e^12*f^6 + 144*a^6*b^12*c^3*e^12*f^6 - 3264*a^7*b^10*c^4*e^12*f^6 + 30112*a^8*b^8
*c^5*e^12*f^6 - 143360*a^9*b^6*c^6*e^12*f^6 + 365568*a^10*b^4*c^7*e^12*f^6 - 458752*a^11*b^2*c^8*e^12*f^6) + 2
04800*a^12*c^9*d*e^11*f^6 + 144*a^6*b^12*c^3*d*e^11*f^6 - 3264*a^7*b^10*c^4*d*e^11*f^6 + 30112*a^8*b^8*c^5*d*e
^11*f^6 - 143360*a^9*b^6*c^6*d*e^11*f^6 + 365568*a^10*b^4*c^7*d*e^11*f^6 - 458752*a^11*b^2*c^8*d*e^11*f^6)*1i)
/((-(9*b^13 - 9*b^4*(-(4*a*c - b^2)^9)^(1/2) + 26880*a^6*b*c^6 + 2077*a^2*b^9*c^2 - 10656*a^3*b^7*c^3 + 30240*
a^4*b^5*c^4 - 44800*a^5*b^3*c^5 - 25*a^2*c^2*(-(4*a*c - b^2)^9)^(1/2) - 213*a*b^11*c + 51*a*b^2*c*(-(4*a*c - b
^2)^9)^(1/2))/(32*(a^5*b^12*e^2*f^4 + 4096*a^11*c^6*e^2*f^4 + 240*a^7*b^8*c^2*e^2*f^4 - 1280*a^8*b^6*c^3*e^2*f
^4 + 3840*a^9*b^4*c^4*e^2*f^4 - 6144*a^10*b^2*c^5*e^2*f^4 - 24*a^6*b^10*c*e^2*f^4)))^(1/2)*((-(9*b^13 - 9*b^4*
(-(4*a*c - b^2)^9)^(1/2) + 26880*a^6*b*c^6 + 2077*a^2*b^9*c^2 - 10656*a^3*b^7*c^3 + 30240*a^4*b^5*c^4 - 44800*
a^5*b^3*c^5 - 25*a^2*c^2*(-(4*a*c - b^2)^9)^(1/2) - 213*a*b^11*c + 51*a*b^2*c*(-(4*a*c - b^2)^9)^(1/2))/(32*(a
^5*b^12*e^2*f^4 + 4096*a^11*c^6*e^2*f^4 + 240*a^7*b^8*c^2*e^2*f^4 - 1280*a^8*b^6*c^3*e^2*f^4 + 3840*a^9*b^4*c^
4*e^2*f^4 - 6144*a^10*b^2*c^5*e^2*f^4 - 24*a^6*b^10*c*e^2*f^4)))^(1/2)*((-(9*b^13 - 9*b^4*(-(4*a*c - b^2)^9)^(
1/2) + 26880*a^6*b*c^6 + 2077*a^2*b^9*c^2 - 10656*a^3*b^7*c^3 + 30240*a^4*b^5*c^4 - 44800*a^5*b^3*c^5 - 25*a^2
*c^2*(-(4*a*c - b^2)^9)^(1/2) - 213*a*b^11*c + 51*a*b^2*c*(-(4*a*c - b^2)^9)^(1/2))/(32*(a^5*b^12*e^2*f^4 + 40
96*a^11*c^6*e^2*f^4 + 240*a^7*b^8*c^2*e^2*f^4 - 1280*a^8*b^6*c^3*e^2*f^4 + 3840*a^9*b^4*c^4*e^2*f^4 - 6144*a^1
0*b^2*c^5*e^2*f^4 - 24*a^6*b^10*c*e^2*f^4)))^(1/2)*(x*(256*a^10*b^13*c^2*e^14*f^10 - 6144*a^11*b^11*c^3*e^14*f
^10 + 61440*a^12*b^9*c^4*e^14*f^10 - 327680*a^13*b^7*c^5*e^14*f^10 + 983040*a^14*b^5*c^6*e^14*f^10 - 1572864*a
^15*b^3*c^7*e^14*f^10 + 1048576*a^16*b*c^8*e^14*f^10) + 1048576*a^16*b*c^8*d*e^13*f^10 + 256*a^10*b^13*c^2*d*e
^13*f^10 - 6144*a^11*b^11*c^3*d*e^13*f^10 + 61440*a^12*b^9*c^4*d*e^13*f^10 - 327680*a^13*b^7*c^5*d*e^13*f^10 +
983040*a^14*b^5*c^6*d*e^13*f^10 - 1572864*a^15*b^3*c^7*d*e^13*f^10) + 192*a^8*b^13*c^2*e^12*f^8 - 4672*a^9*b^
11*c^3*e^12*f^8 + 47360*a^10*b^9*c^4*e^12*f^8 - 256000*a^11*b^7*c^5*e^12*f^8 + 778240*a^12*b^5*c^6*e^12*f^8 -
1261568*a^13*b^3*c^7*e^12*f^8 + 851968*a^14*b*c^8*e^12*f^8) + x*(204800*a^12*c^9*e^12*f^6 + 144*a^6*b^12*c^3*e
^12*f^6 - 3264*a^7*b^10*c^4*e^12*f^6 + 30112*a^8*b^8*c^5*e^12*f^6 - 143360*a^9*b^6*c^6*e^12*f^6 + 365568*a^10*
b^4*c^7*e^12*f^6 - 458752*a^11*b^2*c^8*e^12*f^6) + 204800*a^12*c^9*d*e^11*f^6 + 144*a^6*b^12*c^3*d*e^11*f^6 -
3264*a^7*b^10*c^4*d*e^11*f^6 + 30112*a^8*b^8*c^5*d*e^11*f^6 - 143360*a^9*b^6*c^6*d*e^11*f^6 + 365568*a^10*b^4*
c^7*d*e^11*f^6 - 458752*a^11*b^2*c^8*d*e^11*f^6) - (-(9*b^13 - 9*b^4*(-(4*a*c - b^2)^9)^(1/2) + 26880*a^6*b*c^
6 + 2077*a^2*b^9*c^2 - 10656*a^3*b^7*c^3 + 30240*a^4*b^5*c^4 - 44800*a^5*b^3*c^5 - 25*a^2*c^2*(-(4*a*c - b^2)^
9)^(1/2) - 213*a*b^11*c + 51*a*b^2*c*(-(4*a*c - b^2)^9)^(1/2))/(32*(a^5*b^12*e^2*f^4 + 4096*a^11*c^6*e^2*f^4 +
240*a^7*b^8*c^2*e^2*f^4 - 1280*a^8*b^6*c^3*e^2*f^4 + 3840*a^9*b^4*c^4*e^2*f^4 - 6144*a^10*b^2*c^5*e^2*f^4 - 2
4*a^6*b^10*c*e^2*f^4)))^(1/2)*((-(9*b^13 - 9*b^4*(-(4*a*c - b^2)^9)^(1/2) + 26880*a^6*b*c^6 + 2077*a^2*b^9*c^2
- 10656*a^3*b^7*c^3 + 30240*a^4*b^5*c^4 - 44800*a^5*b^3*c^5 - 25*a^2*c^2*(-(4*a*c - b^2)^9)^(1/2) - 213*a*b^1
1*c + 51*a*b^2*c*(-(4*a*c - b^2)^9)^(1/2))/(32*(a^5*b^12*e^2*f^4 + 4096*a^11*c^6*e^2*f^4 + 240*a^7*b^8*c^2*e^2
*f^4 - 1280*a^8*b^6*c^3*e^2*f^4 + 3840*a^9*b^4*c^4*e^2*f^4 - 6144*a^10*b^2*c^5*e^2*f^4 - 24*a^6*b^10*c*e^2*f^4
)))^(1/2)*((-(9*b^13 - 9*b^4*(-(4*a*c - b^2)^9)^(1/2) + 26880*a^6*b*c^6 + 2077*a^2*b^9*c^2 - 10656*a^3*b^7*c^3
+ 30240*a^4*b^5*c^4 - 44800*a^5*b^3*c^5 - 25*a^2*c^2*(-(4*a*c - b^2)^9)^(1/2) - 213*a*b^11*c + 51*a*b^2*c*(-(
4*a*c - b^2)^9)^(1/2))/(32*(a^5*b^12*e^2*f^4 + 4096*a^11*c^6*e^2*f^4 + 240*a^7*b^8*c^2*e^2*f^4 - 1280*a^8*b^6*
c^3*e^2*f^4 + 3840*a^9*b^4*c^4*e^2*f^4 - 6144*a^10*b^2*c^5*e^2*f^4 - 24*a^6*b^10*c*e^2*f^4)))^(1/2)*(x*(256*a^
10*b^13*c^2*e^14*f^10 - 6144*a^11*b^11*c^3*e^14*f^10 + 61440*a^12*b^9*c^4*e^14*f^10 - 327680*a^13*b^7*c^5*e^14
*f^10 + 983040*a^14*b^5*c^6*e^14*f^10 - 1572864*a^15*b^3*c^7*e^14*f^10 + 1048576*a^16*b*c^8*e^14*f^10) + 10485
76*a^16*b*c^8*d*e^13*f^10 + 256*a^10*b^13*c^2*d*e^13*f^10 - 6144*a^11*b^11*c^3*d*e^13*f^10 + 61440*a^12*b^9*c^
4*d*e^13*f^10 - 327680*a^13*b^7*c^5*d*e^13*f^10 + 983040*a^14*b^5*c^6*d*e^13*f^10 - 1572864*a^15*b^3*c^7*d*e^1
3*f^10) - 192*a^8*b^13*c^2*e^12*f^8 + 4672*a^9*b^11*c^3*e^12*f^8 - 47360*a^10*b^9*c^4*e^12*f^8 + 256000*a^11*b
^7*c^5*e^12*f^8 - 778240*a^12*b^5*c^6*e^12*f^8 + 1261568*a^13*b^3*c^7*e^12*f^8 - 851968*a^14*b*c^8*e^12*f^8) +
x*(204800*a^12*c^9*e^12*f^6 + 144*a^6*b^12*c^3*e^12*f^6 - 3264*a^7*b^10*c^4*e^12*f^6 + 30112*a^8*b^8*c^5*e^12
*f^6 - 143360*a^9*b^6*c^6*e^12*f^6 + 365568*a^10*b^4*c^7*e^12*f^6 - 458752*a^11*b^2*c^8*e^12*f^6) + 204800*a^1
2*c^9*d*e^11*f^6 + 144*a^6*b^12*c^3*d*e^11*f^6 - 3264*a^7*b^10*c^4*d*e^11*f^6 + 30112*a^8*b^8*c^5*d*e^11*f^6 -
143360*a^9*b^6*c^6*d*e^11*f^6 + 365568*a^10*b^4*c^7*d*e^11*f^6 - 458752*a^11*b^2*c^8*d*e^11*f^6) + 128000*a^1
0*c^9*e^10*f^4 + 504*a^6*b^8*c^5*e^10*f^4 - 8112*a^7*b^6*c^6*e^10*f^4 + 48704*a^8*b^4*c^7*e^10*f^4 - 129280*a^
9*b^2*c^8*e^10*f^4))*(-(9*b^13 - 9*b^4*(-(4*a*c - b^2)^9)^(1/2) + 26880*a^6*b*c^6 + 2077*a^2*b^9*c^2 - 10656*a
^3*b^7*c^3 + 30240*a^4*b^5*c^4 - 44800*a^5*b^3*c^5 - 25*a^2*c^2*(-(4*a*c - b^2)^9)^(1/2) - 213*a*b^11*c + 51*a
*b^2*c*(-(4*a*c - b^2)^9)^(1/2))/(32*(a^5*b^12*e^2*f^4 + 4096*a^11*c^6*e^2*f^4 + 240*a^7*b^8*c^2*e^2*f^4 - 128
0*a^8*b^6*c^3*e^2*f^4 + 3840*a^9*b^4*c^4*e^2*f^4 - 6144*a^10*b^2*c^5*e^2*f^4 - 24*a^6*b^10*c*e^2*f^4)))^(1/2)*
2i - atan(((-(9*b^13 + 9*b^4*(-(4*a*c - b^2)^9)^(1/2) + 26880*a^6*b*c^6 + 2077*a^2*b^9*c^2 - 10656*a^3*b^7*c^3
+ 30240*a^4*b^5*c^4 - 44800*a^5*b^3*c^5 + 25*a^2*c^2*(-(4*a*c - b^2)^9)^(1/2) - 213*a*b^11*c - 51*a*b^2*c*(-(
4*a*c - b^2)^9)^(1/2))/(32*(a^5*b^12*e^2*f^4 + 4096*a^11*c^6*e^2*f^4 + 240*a^7*b^8*c^2*e^2*f^4 - 1280*a^8*b^6*
c^3*e^2*f^4 + 3840*a^9*b^4*c^4*e^2*f^4 - 6144*a^10*b^2*c^5*e^2*f^4 - 24*a^6*b^10*c*e^2*f^4)))^(1/2)*((-(9*b^13
+ 9*b^4*(-(4*a*c - b^2)^9)^(1/2) + 26880*a^6*b*c^6 + 2077*a^2*b^9*c^2 - 10656*a^3*b^7*c^3 + 30240*a^4*b^5*c^4
- 44800*a^5*b^3*c^5 + 25*a^2*c^2*(-(4*a*c - b^2)^9)^(1/2) - 213*a*b^11*c - 51*a*b^2*c*(-(4*a*c - b^2)^9)^(1/2
))/(32*(a^5*b^12*e^2*f^4 + 4096*a^11*c^6*e^2*f^4 + 240*a^7*b^8*c^2*e^2*f^4 - 1280*a^8*b^6*c^3*e^2*f^4 + 3840*a
^9*b^4*c^4*e^2*f^4 - 6144*a^10*b^2*c^5*e^2*f^4 - 24*a^6*b^10*c*e^2*f^4)))^(1/2)*((-(9*b^13 + 9*b^4*(-(4*a*c -
b^2)^9)^(1/2) + 26880*a^6*b*c^6 + 2077*a^2*b^9*c^2 - 10656*a^3*b^7*c^3 + 30240*a^4*b^5*c^4 - 44800*a^5*b^3*c^5
+ 25*a^2*c^2*(-(4*a*c - b^2)^9)^(1/2) - 213*a*b^11*c - 51*a*b^2*c*(-(4*a*c - b^2)^9)^(1/2))/(32*(a^5*b^12*e^2
*f^4 + 4096*a^11*c^6*e^2*f^4 + 240*a^7*b^8*c^2*e^2*f^4 - 1280*a^8*b^6*c^3*e^2*f^4 + 3840*a^9*b^4*c^4*e^2*f^4 -
6144*a^10*b^2*c^5*e^2*f^4 - 24*a^6*b^10*c*e^2*f^4)))^(1/2)*(x*(256*a^10*b^13*c^2*e^14*f^10 - 6144*a^11*b^11*c
^3*e^14*f^10 + 61440*a^12*b^9*c^4*e^14*f^10 - 327680*a^13*b^7*c^5*e^14*f^10 + 983040*a^14*b^5*c^6*e^14*f^10 -
1572864*a^15*b^3*c^7*e^14*f^10 + 1048576*a^16*b*c^8*e^14*f^10) + 1048576*a^16*b*c^8*d*e^13*f^10 + 256*a^10*b^1
3*c^2*d*e^13*f^10 - 6144*a^11*b^11*c^3*d*e^13*f^10 + 61440*a^12*b^9*c^4*d*e^13*f^10 - 327680*a^13*b^7*c^5*d*e^
13*f^10 + 983040*a^14*b^5*c^6*d*e^13*f^10 - 1572864*a^15*b^3*c^7*d*e^13*f^10) - 192*a^8*b^13*c^2*e^12*f^8 + 46
72*a^9*b^11*c^3*e^12*f^8 - 47360*a^10*b^9*c^4*e^12*f^8 + 256000*a^11*b^7*c^5*e^12*f^8 - 778240*a^12*b^5*c^6*e^
12*f^8 + 1261568*a^13*b^3*c^7*e^12*f^8 - 851968*a^14*b*c^8*e^12*f^8) + x*(204800*a^12*c^9*e^12*f^6 + 144*a^6*b
^12*c^3*e^12*f^6 - 3264*a^7*b^10*c^4*e^12*f^6 + 30112*a^8*b^8*c^5*e^12*f^6 - 143360*a^9*b^6*c^6*e^12*f^6 + 365
568*a^10*b^4*c^7*e^12*f^6 - 458752*a^11*b^2*c^8*e^12*f^6) + 204800*a^12*c^9*d*e^11*f^6 + 144*a^6*b^12*c^3*d*e^
11*f^6 - 3264*a^7*b^10*c^4*d*e^11*f^6 + 30112*a^8*b^8*c^5*d*e^11*f^6 - 143360*a^9*b^6*c^6*d*e^11*f^6 + 365568*
a^10*b^4*c^7*d*e^11*f^6 - 458752*a^11*b^2*c^8*d*e^11*f^6)*1i + (-(9*b^13 + 9*b^4*(-(4*a*c - b^2)^9)^(1/2) + 26
880*a^6*b*c^6 + 2077*a^2*b^9*c^2 - 10656*a^3*b^7*c^3 + 30240*a^4*b^5*c^4 - 44800*a^5*b^3*c^5 + 25*a^2*c^2*(-(4
*a*c - b^2)^9)^(1/2) - 213*a*b^11*c - 51*a*b^2*c*(-(4*a*c - b^2)^9)^(1/2))/(32*(a^5*b^12*e^2*f^4 + 4096*a^11*c
^6*e^2*f^4 + 240*a^7*b^8*c^2*e^2*f^4 - 1280*a^8*b^6*c^3*e^2*f^4 + 3840*a^9*b^4*c^4*e^2*f^4 - 6144*a^10*b^2*c^5
*e^2*f^4 - 24*a^6*b^10*c*e^2*f^4)))^(1/2)*((-(9*b^13 + 9*b^4*(-(4*a*c - b^2)^9)^(1/2) + 26880*a^6*b*c^6 + 2077
*a^2*b^9*c^2 - 10656*a^3*b^7*c^3 + 30240*a^4*b^5*c^4 - 44800*a^5*b^3*c^5 + 25*a^2*c^2*(-(4*a*c - b^2)^9)^(1/2)
- 213*a*b^11*c - 51*a*b^2*c*(-(4*a*c - b^2)^9)^(1/2))/(32*(a^5*b^12*e^2*f^4 + 4096*a^11*c^6*e^2*f^4 + 240*a^7
*b^8*c^2*e^2*f^4 - 1280*a^8*b^6*c^3*e^2*f^4 + 3840*a^9*b^4*c^4*e^2*f^4 - 6144*a^10*b^2*c^5*e^2*f^4 - 24*a^6*b^
10*c*e^2*f^4)))^(1/2)*((-(9*b^13 + 9*b^4*(-(4*a*c - b^2)^9)^(1/2) + 26880*a^6*b*c^6 + 2077*a^2*b^9*c^2 - 10656
*a^3*b^7*c^3 + 30240*a^4*b^5*c^4 - 44800*a^5*b^3*c^5 + 25*a^2*c^2*(-(4*a*c - b^2)^9)^(1/2) - 213*a*b^11*c - 51
*a*b^2*c*(-(4*a*c - b^2)^9)^(1/2))/(32*(a^5*b^12*e^2*f^4 + 4096*a^11*c^6*e^2*f^4 + 240*a^7*b^8*c^2*e^2*f^4 - 1
280*a^8*b^6*c^3*e^2*f^4 + 3840*a^9*b^4*c^4*e^2*f^4 - 6144*a^10*b^2*c^5*e^2*f^4 - 24*a^6*b^10*c*e^2*f^4)))^(1/2
)*(x*(256*a^10*b^13*c^2*e^14*f^10 - 6144*a^11*b^11*c^3*e^14*f^10 + 61440*a^12*b^9*c^4*e^14*f^10 - 327680*a^13*
b^7*c^5*e^14*f^10 + 983040*a^14*b^5*c^6*e^14*f^10 - 1572864*a^15*b^3*c^7*e^14*f^10 + 1048576*a^16*b*c^8*e^14*f
^10) + 1048576*a^16*b*c^8*d*e^13*f^10 + 256*a^10*b^13*c^2*d*e^13*f^10 - 6144*a^11*b^11*c^3*d*e^13*f^10 + 61440
*a^12*b^9*c^4*d*e^13*f^10 - 327680*a^13*b^7*c^5*d*e^13*f^10 + 983040*a^14*b^5*c^6*d*e^13*f^10 - 1572864*a^15*b
^3*c^7*d*e^13*f^10) + 192*a^8*b^13*c^2*e^12*f^8 - 4672*a^9*b^11*c^3*e^12*f^8 + 47360*a^10*b^9*c^4*e^12*f^8 - 2
56000*a^11*b^7*c^5*e^12*f^8 + 778240*a^12*b^5*c^6*e^12*f^8 - 1261568*a^13*b^3*c^7*e^12*f^8 + 851968*a^14*b*c^8
*e^12*f^8) + x*(204800*a^12*c^9*e^12*f^6 + 144*a^6*b^12*c^3*e^12*f^6 - 3264*a^7*b^10*c^4*e^12*f^6 + 30112*a^8*
b^8*c^5*e^12*f^6 - 143360*a^9*b^6*c^6*e^12*f^6 + 365568*a^10*b^4*c^7*e^12*f^6 - 458752*a^11*b^2*c^8*e^12*f^6)
+ 204800*a^12*c^9*d*e^11*f^6 + 144*a^6*b^12*c^3*d*e^11*f^6 - 3264*a^7*b^10*c^4*d*e^11*f^6 + 30112*a^8*b^8*c^5*
d*e^11*f^6 - 143360*a^9*b^6*c^6*d*e^11*f^6 + 365568*a^10*b^4*c^7*d*e^11*f^6 - 458752*a^11*b^2*c^8*d*e^11*f^6)*
1i)/((-(9*b^13 + 9*b^4*(-(4*a*c - b^2)^9)^(1/2) + 26880*a^6*b*c^6 + 2077*a^2*b^9*c^2 - 10656*a^3*b^7*c^3 + 302
40*a^4*b^5*c^4 - 44800*a^5*b^3*c^5 + 25*a^2*c^2*(-(4*a*c - b^2)^9)^(1/2) - 213*a*b^11*c - 51*a*b^2*c*(-(4*a*c
- b^2)^9)^(1/2))/(32*(a^5*b^12*e^2*f^4 + 4096*a^11*c^6*e^2*f^4 + 240*a^7*b^8*c^2*e^2*f^4 - 1280*a^8*b^6*c^3*e^
2*f^4 + 3840*a^9*b^4*c^4*e^2*f^4 - 6144*a^10*b^2*c^5*e^2*f^4 - 24*a^6*b^10*c*e^2*f^4)))^(1/2)*((-(9*b^13 + 9*b
^4*(-(4*a*c - b^2)^9)^(1/2) + 26880*a^6*b*c^6 + 2077*a^2*b^9*c^2 - 10656*a^3*b^7*c^3 + 30240*a^4*b^5*c^4 - 448
00*a^5*b^3*c^5 + 25*a^2*c^2*(-(4*a*c - b^2)^9)^(1/2) - 213*a*b^11*c - 51*a*b^2*c*(-(4*a*c - b^2)^9)^(1/2))/(32
*(a^5*b^12*e^2*f^4 + 4096*a^11*c^6*e^2*f^4 + 240*a^7*b^8*c^2*e^2*f^4 - 1280*a^8*b^6*c^3*e^2*f^4 + 3840*a^9*b^4
*c^4*e^2*f^4 - 6144*a^10*b^2*c^5*e^2*f^4 - 24*a^6*b^10*c*e^2*f^4)))^(1/2)*((-(9*b^13 + 9*b^4*(-(4*a*c - b^2)^9
)^(1/2) + 26880*a^6*b*c^6 + 2077*a^2*b^9*c^2 - 10656*a^3*b^7*c^3 + 30240*a^4*b^5*c^4 - 44800*a^5*b^3*c^5 + 25*
a^2*c^2*(-(4*a*c - b^2)^9)^(1/2) - 213*a*b^11*c - 51*a*b^2*c*(-(4*a*c - b^2)^9)^(1/2))/(32*(a^5*b^12*e^2*f^4 +
4096*a^11*c^6*e^2*f^4 + 240*a^7*b^8*c^2*e^2*f^4 - 1280*a^8*b^6*c^3*e^2*f^4 + 3840*a^9*b^4*c^4*e^2*f^4 - 6144*
a^10*b^2*c^5*e^2*f^4 - 24*a^6*b^10*c*e^2*f^4)))^(1/2)*(x*(256*a^10*b^13*c^2*e^14*f^10 - 6144*a^11*b^11*c^3*e^1
4*f^10 + 61440*a^12*b^9*c^4*e^14*f^10 - 327680*a^13*b^7*c^5*e^14*f^10 + 983040*a^14*b^5*c^6*e^14*f^10 - 157286
4*a^15*b^3*c^7*e^14*f^10 + 1048576*a^16*b*c^8*e^14*f^10) + 1048576*a^16*b*c^8*d*e^13*f^10 + 256*a^10*b^13*c^2*
d*e^13*f^10 - 6144*a^11*b^11*c^3*d*e^13*f^10 + 61440*a^12*b^9*c^4*d*e^13*f^10 - 327680*a^13*b^7*c^5*d*e^13*f^1
0 + 983040*a^14*b^5*c^6*d*e^13*f^10 - 1572864*a^15*b^3*c^7*d*e^13*f^10) + 192*a^8*b^13*c^2*e^12*f^8 - 4672*a^9
*b^11*c^3*e^12*f^8 + 47360*a^10*b^9*c^4*e^12*f^8 - 256000*a^11*b^7*c^5*e^12*f^8 + 778240*a^12*b^5*c^6*e^12*f^8
- 1261568*a^13*b^3*c^7*e^12*f^8 + 851968*a^14*b*c^8*e^12*f^8) + x*(204800*a^12*c^9*e^12*f^6 + 144*a^6*b^12*c^
3*e^12*f^6 - 3264*a^7*b^10*c^4*e^12*f^6 + 30112*a^8*b^8*c^5*e^12*f^6 - 143360*a^9*b^6*c^6*e^12*f^6 + 365568*a^
10*b^4*c^7*e^12*f^6 - 458752*a^11*b^2*c^8*e^12*f^6) + 204800*a^12*c^9*d*e^11*f^6 + 144*a^6*b^12*c^3*d*e^11*f^6
- 3264*a^7*b^10*c^4*d*e^11*f^6 + 30112*a^8*b^8*c^5*d*e^11*f^6 - 143360*a^9*b^6*c^6*d*e^11*f^6 + 365568*a^10*b
^4*c^7*d*e^11*f^6 - 458752*a^11*b^2*c^8*d*e^11*f^6) - (-(9*b^13 + 9*b^4*(-(4*a*c - b^2)^9)^(1/2) + 26880*a^6*b
*c^6 + 2077*a^2*b^9*c^2 - 10656*a^3*b^7*c^3 + 30240*a^4*b^5*c^4 - 44800*a^5*b^3*c^5 + 25*a^2*c^2*(-(4*a*c - b^
2)^9)^(1/2) - 213*a*b^11*c - 51*a*b^2*c*(-(4*a*c - b^2)^9)^(1/2))/(32*(a^5*b^12*e^2*f^4 + 4096*a^11*c^6*e^2*f^
4 + 240*a^7*b^8*c^2*e^2*f^4 - 1280*a^8*b^6*c^3*e^2*f^4 + 3840*a^9*b^4*c^4*e^2*f^4 - 6144*a^10*b^2*c^5*e^2*f^4
- 24*a^6*b^10*c*e^2*f^4)))^(1/2)*((-(9*b^13 + 9*b^4*(-(4*a*c - b^2)^9)^(1/2) + 26880*a^6*b*c^6 + 2077*a^2*b^9*
c^2 - 10656*a^3*b^7*c^3 + 30240*a^4*b^5*c^4 - 44800*a^5*b^3*c^5 + 25*a^2*c^2*(-(4*a*c - b^2)^9)^(1/2) - 213*a*
b^11*c - 51*a*b^2*c*(-(4*a*c - b^2)^9)^(1/2))/(32*(a^5*b^12*e^2*f^4 + 4096*a^11*c^6*e^2*f^4 + 240*a^7*b^8*c^2*
e^2*f^4 - 1280*a^8*b^6*c^3*e^2*f^4 + 3840*a^9*b^4*c^4*e^2*f^4 - 6144*a^10*b^2*c^5*e^2*f^4 - 24*a^6*b^10*c*e^2*
f^4)))^(1/2)*((-(9*b^13 + 9*b^4*(-(4*a*c - b^2)^9)^(1/2) + 26880*a^6*b*c^6 + 2077*a^2*b^9*c^2 - 10656*a^3*b^7*
c^3 + 30240*a^4*b^5*c^4 - 44800*a^5*b^3*c^5 + 25*a^2*c^2*(-(4*a*c - b^2)^9)^(1/2) - 213*a*b^11*c - 51*a*b^2*c*
(-(4*a*c - b^2)^9)^(1/2))/(32*(a^5*b^12*e^2*f^4 + 4096*a^11*c^6*e^2*f^4 + 240*a^7*b^8*c^2*e^2*f^4 - 1280*a^8*b
^6*c^3*e^2*f^4 + 3840*a^9*b^4*c^4*e^2*f^4 - 6144*a^10*b^2*c^5*e^2*f^4 - 24*a^6*b^10*c*e^2*f^4)))^(1/2)*(x*(256
*a^10*b^13*c^2*e^14*f^10 - 6144*a^11*b^11*c^3*e^14*f^10 + 61440*a^12*b^9*c^4*e^14*f^10 - 327680*a^13*b^7*c^5*e
^14*f^10 + 983040*a^14*b^5*c^6*e^14*f^10 - 1572864*a^15*b^3*c^7*e^14*f^10 + 1048576*a^16*b*c^8*e^14*f^10) + 10
48576*a^16*b*c^8*d*e^13*f^10 + 256*a^10*b^13*c^2*d*e^13*f^10 - 6144*a^11*b^11*c^3*d*e^13*f^10 + 61440*a^12*b^9
*c^4*d*e^13*f^10 - 327680*a^13*b^7*c^5*d*e^13*f^10 + 983040*a^14*b^5*c^6*d*e^13*f^10 - 1572864*a^15*b^3*c^7*d*
e^13*f^10) - 192*a^8*b^13*c^2*e^12*f^8 + 4672*a^9*b^11*c^3*e^12*f^8 - 47360*a^10*b^9*c^4*e^12*f^8 + 256000*a^1
1*b^7*c^5*e^12*f^8 - 778240*a^12*b^5*c^6*e^12*f^8 + 1261568*a^13*b^3*c^7*e^12*f^8 - 851968*a^14*b*c^8*e^12*f^8
) + x*(204800*a^12*c^9*e^12*f^6 + 144*a^6*b^12*c^3*e^12*f^6 - 3264*a^7*b^10*c^4*e^12*f^6 + 30112*a^8*b^8*c^5*e
^12*f^6 - 143360*a^9*b^6*c^6*e^12*f^6 + 365568*a^10*b^4*c^7*e^12*f^6 - 458752*a^11*b^2*c^8*e^12*f^6) + 204800*
a^12*c^9*d*e^11*f^6 + 144*a^6*b^12*c^3*d*e^11*f^6 - 3264*a^7*b^10*c^4*d*e^11*f^6 + 30112*a^8*b^8*c^5*d*e^11*f^
6 - 143360*a^9*b^6*c^6*d*e^11*f^6 + 365568*a^10*b^4*c^7*d*e^11*f^6 - 458752*a^11*b^2*c^8*d*e^11*f^6) + 128000*
a^10*c^9*e^10*f^4 + 504*a^6*b^8*c^5*e^10*f^4 - 8112*a^7*b^6*c^6*e^10*f^4 + 48704*a^8*b^4*c^7*e^10*f^4 - 129280
*a^9*b^2*c^8*e^10*f^4))*(-(9*b^13 + 9*b^4*(-(4*a*c - b^2)^9)^(1/2) + 26880*a^6*b*c^6 + 2077*a^2*b^9*c^2 - 1065
6*a^3*b^7*c^3 + 30240*a^4*b^5*c^4 - 44800*a^5*b^3*c^5 + 25*a^2*c^2*(-(4*a*c - b^2)^9)^(1/2) - 213*a*b^11*c - 5
1*a*b^2*c*(-(4*a*c - b^2)^9)^(1/2))/(32*(a^5*b^12*e^2*f^4 + 4096*a^11*c^6*e^2*f^4 + 240*a^7*b^8*c^2*e^2*f^4 -
1280*a^8*b^6*c^3*e^2*f^4 + 3840*a^9*b^4*c^4*e^2*f^4 - 6144*a^10*b^2*c^5*e^2*f^4 - 24*a^6*b^10*c*e^2*f^4)))^(1/
2)*2i - ((x*(3*b^3*d - 20*a*c^2*d^3 + 6*b^2*c*d^3 - 11*a*b*c*d))/(a*(a*b^2 - 4*a^2*c)) - (x^4*(10*a*c^2*e^3 -
3*b^2*c*e^3))/(2*a*(a*b^2 - 4*a^2*c)) - (2*x^3*(10*a*c^2*d*e^2 - 3*b^2*c*d*e^2))/(a*(a*b^2 - 4*a^2*c)) + (2*a*
b^2 - 8*a^2*c + 3*b^3*d^2 - 10*a*c^2*d^4 + 3*b^2*c*d^4 - 11*a*b*c*d^2)/(2*a*e*(a*b^2 - 4*a^2*c)) + (x^2*(3*b^3
*e - 60*a*c^2*d^2*e + 18*b^2*c*d^2*e - 11*a*b*c*e))/(2*a*(a*b^2 - 4*a^2*c)))/(x^2*(10*c*d^3*e^2*f^2 + 3*b*d*e^
2*f^2) + x*(a*e*f^2 + 3*b*d^2*e*f^2 + 5*c*d^4*e*f^2) + x^3*(b*e^3*f^2 + 10*c*d^2*e^3*f^2) + b*d^3*f^2 + c*d^5*
f^2 + a*d*f^2 + c*e^5*f^2*x^5 + 5*c*d*e^4*f^2*x^4)